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ABSTRACT 

Gamma ray bursts can potentially be used as distance indicators, providing the pos- 
sibility of extending the Hubble diagram to redshifts ^ 7. Here we follow the analysis 
of Schaefer (2007), with the aim of distinguishing the timescape cosmological model 
from the ACDM model by means of the additional leverage provided by GRBs in the 
range 2 < z < 7. We find that the timescape model fits the GRB sample slightly 
better than the ACDM model, but that the systematic uncertainties are still too little 
understood to distinguish the models. 
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1 INTRODUCTION 

The timescape (TS) model is an inhomogeneous cosmolog- 
ical model that explains the apparently accelerated cosmic 
expansion first observed in the supernova luminosity dis- 
tances as an artifact of gradients in gravitational energy 
between gravitationally bound systems and the interven- 
ing negatively curved voids. In the dynamical spacetime of 
general relativity, this leads to a variance in the calibra- 
tion of the clock rates of ideal observers who fit average 
smoothed-out geometries to the underlying inhomogeneous 
matter distribution. The TS model agrees closely with the 
ACDM model over the range of scales probed by the su- 
pernova data (Leith et al. 2008), with certain qualifications: 
parameter values obtained by minimizing y^ fits to the TS 
Hubble curve depend significantly on the process used to 
reduce the SN la light curves (Smale & Wiltshire 2011). 

The current state of knowledge of systematic uncertain- 
ties in the SN la data precludes discrimination between the 
TS and ACDM models using SNe la (Smale & Wiltshire 
2011). In fact, calculation of the eflfective comoving distance 
HoD{z) shows that in the redshift range probed by SNe la 
there is little to distinguish between the TS model with the 
best-fit value for the present void fraction /vo = 0.762 from 
the Gold dataset of Riess et al. (2007). Wihshire (2009) has 
noted that over different redshift ranges HoD{z) for the TS 
model closely approximates HqD{z) for spatially flat ACDM 
models with different values of fimo and Hao • It is thus seen 
to interpolate between different ACDM models as the red- 
shift is varied (see Fig. 1). Fig. 1 shows that between z ~ 2 
and 2; ~ 6, the TS HoD{z) crosses from coinciding closely 
with the best fit line from the SNe la only to that predicted 
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by the best fit to WMAP, BAG and the SNe la. In princi- 
ple. Gamma Ray Bursts (GRBs), which probe this redshift 
range, could distinguish the TS and ACDM models in this 
redshift range, although their use as distance indicators is 
far from established. 

This paper will establish that the TS model is also sup- 
ported by the current GRB data (Schaefer 2007), but that, 
as one might expect from the SNe la results, the uncertain- 
ties in the data are as yet too large to distinguish the models 
in the redshift range 2 < z < 6. 

The paper is organized as follows. Section 2 explains 
the method of "standardizing" the GRBs for their use as 
distance indicators, and gives a brief derivation of the TS 
luminosity distance. Section 3 describes the results, before 
a discussion and conclusion are presented in Section 4. 



2 BACKGROUND 

2.1 The timescape model 

In keeping with current observations that the large-scale cos- 
mic structure consists of voids of average diameter ~ 30 /i~^ 
Mpc (Hoyle & Vogeley 2004; Pan et al. 2011) separated 
and threaded by walls and filaments containing clusters 
of galaxies, the timescape model is based on a differentia- 
tion of the Universe into gravitationally bound spatially flat 
wall regions and negatively curved voids. There is a con- 
sequent small backreaction (^ 5% as a normalized energy 
density) which nevertheless leads to significant cosmologi- 
cal effects over cosmological timescales (Wiltshire 2007a). 
At late epochs the construction of a single smoothed-out 
geometry becomes problematic when the underlying geom- 
etry varies. The TS model is based on the assumption that 
different equivalent descriptions of a smoothed-out average 



Smale 



Figure 1. Effective comoving distance as a function of redshift 
for various spatially flat models (dotted lines) and for the TS 
model with /vo = 0.762 (solid line). Parameter values for the 
dotted Unes are (i) f^ni = 0.249 (best-fit to WMAP only); (ii) 
Qm = 0.279 (joint best-fit to WMAP, BAO and SNe la); (iii) 
Qm = 0.34 (best-fit to Riess et al. (2007) SNe la only). After 
Wiltshire (2009). 
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geometry can be given, but these descriptions will vary be- 
tween canonical observers who each assume that the average 
geometry has the same spatial curvature as the locally de- 
termined geometry. Differences in the calibration of rulers 
and clocks grow cumulatively as the variance in spatial ge- 
ometry grows, and these must be taken into account when 
reconstructing the expansion history of the universe from 
information on null geodesies. 

As observers in galaxies, our local average geometry, 
assumed to be spatially flat with scale factor aw, is given by 



ds 



:yj = -dr^ + aw(r)[dj7w -|- jy^dfi^]. 



(1) 



Finite infinity (Ellis 1984), denoted by fi, demarcates the 
boundary between gravitationally bound and unbound sys- 
tems (Wiltshire 2007a). A similar expression defines the neg- 
atively curved geometry at the centre of a void, with an 
appropriate void time parameter Tv and scale factor Ov 
The volume averaged scale factor 



/viflv + (1 - /v,)a„ 



(2) 



where /v^ ^ 1 is the initial void fraction, evolves according 
to an averaging of the Einstein equations for an inhomo- 
geneous dust cosmology due to Buchert (2000). The corre- 
sponding averaged geometry, in terms of the proper time t 
at a volume average position in freely expanding space, has 
the form 



ds" 



-dr + a{t)Af) + A{f), i)df7^ 



(3) 



where the area function A is defined by an average over the 
particle horizon volume (Wiltshire 2007a) . The local time in 
a wall is related to t by the phenomenological lapse param- 

eter7=|i. 



A single geometry that relates metrics (1) and (3) is con- 
structed by matching the radial null geodesies of wall and 
volume average geometries sharing a common centre (Wilt- 
shire (2007a), § 5.2). Along the radial null geodesies, the line 
elements of the two geometries are simply related by a con- 
formal factor. Once the metric is extended to cosmological 
scales, instead of (3) wall observers describe the large-scale 
universe by the effective metric 



ds" = -dr" + ^\Aff + ri{fj, t)dn% 



(4) 



where r™ = 7(1 - /v)^/-^(l - /vj"^/^r?w(f;, t). Metric (4) 
is the "dressed" geometry that arises when we attempt to 
reconcile our position as observer within a finite infinity re- 
gion with observations of objects at cosmological distances. 
In general, observers (localised within finite infinity regions) 
cannot assume that their measurements of cosmological pa- 
rameters correspond to the global average values. However, 
due to the existence of a scale of statistical homogeneity, the 
values of locally measured cosmological parameters should 
converge towards their global average values as the averag- 
ing volume increases. 

Volume average ("bare") parameters, referred to met- 
ric (3), differ from dressed parameters. The dressed matter 
density is flm ~ 7"^nm0, and the Hubble parameter of met- 
ric (4) is related to the bare Hubble parameter according 
to 



H = ^H- —7. 
' dt ' 



(5) 



It is a feature of the TS model that the variance of pa- 
rameters calculated from observations of nearby objects (on 
scales < 100 h~^ Mpc) will be relatively large because such 
objects lie within the scale of statistical homogeneity. Since 
the volume of space is dominated by voids — of typical di- 
ameter ~ 30 h^^ Mpc — which appear to expand faster 
than the walls, observers in a typical galaxy looking on a 
typical line of sight through local voids will infer that their 
local universe is expanding faster than the global average. 
Eventually a typical line of sight will intersect a sufficient 
number of walls as well as voids to approach the global av- 
erage, which does not change by sampling on ever larger 
scales. The transition scale from from large to small vari- 
ance in the expansion must be larger than the diameter of 
the dominant voids, and is referred to as the scale of sta- 
tistical homogeneity. It is expected to be comparable to the 
BAO scale, ~ 100 h'^ Mpc (Wiltshire 2008). Thus any ob- 
server in a galaxy will typically see a "Hubble bubble" on 
scales < 100 h~^ Mpc. 

In terms of volume average time, t, the luminosity dis- 
tance for wall observers is (Wiltshire 2007b) 



dL = ao{l + z)r^ 



(6) 



The Buchert equations have an exact general solution which 
admits a particular late-time attractor solution, to which 
the general solution converges to within 1% by redshift 
z ~ 37 (Wiltshire 2007b). For this tracker solution, eq. (6) 
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becomes (Wiltshire 2009) 



HodL = {l + zf{Hotf^^ f 
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where y^ = Hot and &» = 2(1 - /vo)(2 + /vo)/(9/vo). In 
Fig. 1, D = dL/{l + z), and the dressed Hubble constant 
Ho, is related to the bare Hubble constant Hq by 



Ho = (4/vo + /vo + 4)/[2(2 + f.o)]Ho. 



(8) 



The dressed Hubble constant is the one whose value 
should coincide with that which is conventionally deter- 
mined on scales greater than the scale of statistical homo- 
geneity, z > 0.033. Using the "Gold" SN la dataset of Riess 
et al. (2007), Leith et al. (2008) found the dressed Hubble 
constant to be Ho = 61. 7^13 km sec^^Mpc^^, (and the 
bare Hubble constant Hq = 48.2 ± 2.6km sec~^ Mpc~^). 
However, since the supernova data magnitudes depend on 
an overall normalization determined from the local distance 
ladder, they cannot be used to determine the Hubble con- 
stant alone. Joint estimates of Hg and /vo which fit both 
the angular diameter distance of the sound horizon in the 
CMB anisotropy data and the baryon acoustic oscillation 
scale in galaxy clustering statistics do provide an inde- 
pendent constraint on the value of H^, however, and on 
this basis Leith et al. (2008) find the range of values of 
the dressed Hubble constant to be roughly constrained to 
lie in the interval 57 < -ff„ < 68 km sec~^Mpc^^. This is 
lower than the S/fgES best estimate of Riess et al. (2011) 
but that survey relies on the calibration of the distance 
ladder using objects that lie within the scale of statisti- 
cal homogeneity, and this may involve complicated system- 
atics given that higher values of Hg are expected below 
the statistical homogeneity scale. Estimates of H^ which 
do not rely on calibration with nearby objects are often 
somewhat lower. For example, Courbin et al. (2010) find 
Ho = 62^4 km sec^^ Mpc^^ using the time delay from 
strong gravitational lensing of quasars, and Beutler et al. 
(2011) estimate Ho = 67 ±3.2 km sec"^ Mpc"^ using the 
WMAP sound horizon-calibrated BAO signal in the 6dF 
galaxy survey. These measurements indicate that a value of 
Ho consistent with the TS model is still to be obtained once 
systematic errors on distance determinations are reduced. 

With the tracker solution, the bare densities flmo and 
Oko can be written in terms of the present void fraction. 
In particular, the dressed matter density, measured by wall 
observers, for which the numerical value is most likely to be 
similar to that of a FLRW model, will be 



fin 



-(2 + /vo) fimO 
O 



-(l-/vo)(2 + /v(,) 



(9) 



for the tracker solution. 

Since there is no nearby GRB sample, there is no GRB 
calibration of Hq, and we work with relative distances only. 
The fitting process therefore results in a best-fit value for 
the single parameter /vo. 



2.2 GRB data reduction method 

In this paper I will use the sample of 69 GRBs selected 
by Schaefer (2007) (henceforward S07) as having sufficient 
light curve data to compute their placement on a Hubble 
diagram. 

GRBs are not standard candles, since their luminosi- 
ties span several orders of magnitude (whether one assumes 
coUimated or isotropic emission). However, there are ongo- 
ing attempts to "standardize" GRBs given their promise for 
cosmology: they occur at higher redshifts than any estab- 
lished standard candles, and radiation in the gamma band 
(^10 keV) is not subject to the same limitations due to 
dust extinction as the optical band (Ghirlanda, Ghisellini 
& Firmani 2006). Certain GRB light curve parameters have 
been found to correlate with each other, ofi'ering the pos- 
sibility of computing a magnitude, much in the same way 
as the Phillips stretch-luminosity relation is used to reduce 
scatter in the SN la Hubble diagram. 

Schaefer (2007) uses four light curve parameters that 
correlate with the luminosity: (1) the lag time riag is the 
time shift between the hard and the soft light curves; (2) 
the light curve variability V is the normalized variance of the 
light curve around a smoothed version of that light curve; 
(3) the peak energy Epea.k is the photon energy at which 
the vF,^ spectrum is brightest; and (4) the minimum rise 
time trt is the shortest time over which the light curve 
rises by half the peak fiux of the pulse. A fifth correlation, 
Speak — E^, relates the peak energy of the light curve to to- 
tal photon energy emitted by the burst. This is the tightest 
of the correlations (Ghirlanda et al. 2004), but it requires 
measurement of a jet break time by which the measured 
(isotropic-equivalent) energy can be corrected for the coUi- 
mation. Along with these luminosity indicators a peak fiux 
P is also measured for a wide range of bandpasses. A bolo- 
metric fiux (or fluence) can be calculated by extrapolating 
to high and low energies using the well-known broken power 
law of Band et al. (1993) for the GRB spectrum and in- 
tegrating over all energies. This brings consistency to the 
brightnesses, and given a cosmological model, permits the 
calculation of an isotropic luminosity 



L = 47rdj^Pboi< 



(10) 



The algorithm goes as follows. The luminosity indicator 
is the independent variable, and from eq. (10) is obtained the 
Y-coordinate. A linear fit to the logarithms of these quanti- 
ties gives an empirical relationship between the luminosity 
indicator and the luminosity. For this fit, we use the bisec- 
tor of the two ordinary least squares fits: that of X against 
Y, and then vice versa (Isobe et al. 1990). We can then use 
this relationship to calculate a theoretical luminosity curve 
for each indicator, based on the luminosity distance obtained 
from eq. (10). In cases where a jet break has been measured, 
Speak is related to the coUimation-corrected energy E^, by 



E^ = 47rdiS'boio(l - cos6'jet)(l + zY 



(11) 



for jet opening angle 6jct and bolometric fluence S'boio- The 
uncertainties in the F-axis quantities log L and log E^ are 
obtained from the uncertainties in the X-axis quantities in 
the standard way. Because the physics of the GRB explo- 
sions is not completely understood, the correlations contain 
some scatter over and above the measurement noise. To ac- 
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count for this, an additional intrinsic uncertainty is esti- 
mated such that the reduced x^ of the indicator-luminosity 
calibration curve is unity. The best-fit lines for these rela- 
tions are given along with their uncertainties in Appendix A. 
From each calculated L or E^, we then recalculate a 
luminosity distance via (10) or (11), from which we obtain 
a distance modulus in the standard way: fi — 5 log dL — 25 
for dh in Mpc. The propagated uncertainties are (Schaefer 
2007) 

1.086crp,„ 



o"m = (2.5criogL)^ + 



(12) 



-Pbolo 

if the bolometric flux Pboio is used, or, if the bolometric 
fluence is used, 



(J^ = (2.5criog£;.^) 



1.086as^ 

•^bolo 



+ 



1.086aFb,„ 

-^bcam 



i^jct) is calculated 



-) , (13) 



where the beam factor Fbcam = (1 " 
from the jet break time tjot. 

Finally, we take a weighted average of all the five dif- 
ferent distance moduli: 



^ = 






with the uncertainty 



(14) 



(15) 



We avoid circularity by performing a simultaneous fit of 
both the cosmology and the luminosity relations (Ghirlanda 
et al. 2004; Schaefer 2007) — i.e. the luminosity relations are 
part of the model. The value of Ho here is arbitrary, since its 
variation changes the Hubble line and the luminosity cali- 
bration of the data in an identical way, resulting merely in a 
change in the overall normalization of the Hubble diagram. 
GRBs do not occur in the local universe, so calibrating the 
GRB Hubble diagram to a value of Hq with any accuracy 
is not possible. This is different to the SN la case, in which 
the calibration of light curve and stretch can be done model- 
independently with nearby SNe la, and then extrapolated to 
objects at higher redshifts. However, regardless of the nor- 
malization, the shape of the curve in the Hubble diagram 
depends solely on /vo- This means that for a range of values 
of /vo, here between 0.0 and 1.0, we calibrate the luminos- 
ity relation and compute the placement of the GRBs on the 
Hubble diagram, and calculate a corresponding range of x^ 
values. The favoured value for /vo is that for which the x^ 
is minimized. 



3 RESULTS 

We show the results of the linear regressions in figures 2- 
6 in black. The intercept a and slope b of the TS model 
calibration line are shown in each figure. For comparison, the 
ACDM-calibrated data points (for a flat Friedmann model 
with Qm ~ 0.27 and w — —1, as calculated in S07) are shown 
in grey. 

The timescape model produces calibrations that are 
within la of the ACDM model in each case. In fact, the TS 
model regression parameters match those of the concordance 
model more closely than regression parameters calculated 
from the variable dark energy equation of state cosmology 



Figure 2. Power law relation between lag time Ti^g, corrected to 
the GRB rest frame, and isotropic luminosity, for 38 GRBs. The 
Icr measurement uncertainties are used for the error bars. The 
ACDM fit as calculated in Schaefer (2007) is shown in grey. The 
intercept a and the slope b for the TS calibration arc shown on the 
plot, and the equation of the best-fit line and the expected uncer- 
tainty in the luminosity so calculated is given in Appendix A(i). 
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Figure 3. Bisector fit of the Variability-Luminosity relation for 
51 GRBs. Larger measurement uncertainties in this relation mean 
it carries less weight in the final luminosity average. The TS inter- 
cept a and the slope b are shown, and the equation of the best-fit 
line and the expected uncertainty in the luminosity so calculated 
are given in Appendix A(ii). 




log[V(1+z)] 

of Riess et al. (2004) (wo = -1.31, w' = 1.48), computed in 
S07 to assess the dependence of the calibration on the input 
cosmology. 

The resulting Hubble diagram for the TS model is 
shown in fig. 7. In the ACDM case, with the "concordance" 
value Qrn ~ 0.27, we obtain a reduced x^ of 1-05 as in 
S07. The parameter values that minimize the HD x^ ^^^ 
Qm = 0.2llo^^ for the ACDM model\ shown in grey, and 



1 This coincides within a standard deviation with Q„ 
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Figure 4. Bisector fit of GRB isotropic luminosity to the peak 
energy -Epcaki corrected to the rest frame of the GRB. A'^ = 64. 
The TS intercept a and the slope b, and the equation of the best-fit 
line and the expected uncertainty in the luminosity so calculated 
arc given in Appendix A(iii). 



Figure 6. Minimum rise time-Luminosity relation for 62 GRBs. 
The TS intercept a and the slope b are shown, and the equation 
of the best-fit line and the expected uncertainty in the luminosity 
so calculated are given in Appendix A(v). 
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Figure 5. -Epcak ~ E-y relation for 27 GRBs. This is the tightest 
of the five power relations but there are fewer data points, since 
calculation of E^ requires identification and measurement of a jet 
break. The TS intercept a and the slope b are shown on the plot, 
and the equation of the best-fit line and the expected uncertainty 
in the E^ so calculated are given in Appendix A(iv). 
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/vo ~ 0.84l(5'2i for the timescape model, for which the re- 
duced x^ was 1.04 for 68 dof, shown in black. This present 
void fraction corresponds to a matter density as measured 
by wall observers via eq. (9) of f^m = 0.231q'2o- However, 
note that there is no a priori reason why the ACDM and 
TS values for f2m should coincide, since the role of this pa- 
rameter in each theory is different. 

The TS model fits the GRB Hubble diagram slightly 
better (lower x^) than the ACDM model. The correspond- 



Figure 7. Hubble diagram for the 69 GRBs of S07. The ACDM 
diagram is shown in grey, and the TS diagram is shown in black. 
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found in S07, which was found by marginalizing over 



0.39; 

the slopes and intercepts of the luminosity relations. 



Redshift (z) 



ing Bayes factor InS = 0.18 indicates Bayesian evidence 
in favour of the timescape model that is "not worth more 
than a bare mention" according to the Jeffreys scale (Kass 
& Raftery 1995). This is apparent, since the competing pre- 
dictions of the models lie well within the range spanned by 
the measurement errors, let alone the systematics, so it can 
only be concluded that GRB cosmology is not yet precise 
enough to distinguish between these models. 

By contrast, preliminary investigations by Schaefer 
(2008) indicate that certain modified gravity models and 
particular exotic forms of dark energy (the Chaplygin gas) 
provide much poorer fits to the GRB data than the stan- 
dard ACDM model. Amongst the alternatives to the ACDM 
model, the TS model therefore enjoys a degree of phe- 
nomenological success which is hard to replicate in a number 
of other scenarios. 
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4 CONCLUSION 

Some issues with the use of GRBs as standard candles 
for constraining cosmological parameters are discussed by 
Ghirlanda, Ghisellini & Firmani (2006), Ghirlanda (2009), 
and Petrosian et al. (2009). In particular, the correlations 
between the isotropic luminosities and the luminosity indica- 
tors are weak in a x sense — physical factors unaccounted for 
are producing large scatter. Strictly speaking, the Epcan^E-, 
correlation is the only relationship with a sufficiently low re- 
duced x^ to admit cosmological parameter estimation, albeit 
with the caveat that the measurement of the jet break time 
assumes a particular fireball model (Ghirlanda 2009). Pet- 
rosian et al. (2009) point out that the luminosity correlations 
are statistical in nature, rather than being, as they should 
ideally, one-to-one relations between uncorrelated quanti- 
ties. This meant, for example, that the comparatively tight 
Speak — E^ correlation found by Ghirlanda et al. (2004) ac- 
tually weakened with the introduction of more data points. 

Systematic uncertainties such as dust extinction and 
evolution constitute considerable limitations to cosmologi- 
cal parameter estimation with SNe la. Many of these un- 
certainties, for example Malmquist bias and gravitational 
lensing, are considered negligible in the redshift range over 
which SNe la occur. For the redshift range over which GRBs 
occur, one would expect that Malmquist bias and lens- 
ing might cause at least some of the scatter in the GRB 
Hubble diagram, but these biases are shown in S07 to be 
negligibly small. Obscuration by dust is not an issue for 
GRBs (Ghirlanda et al. 2004), but selection and evolution ef- 
fects can potentially influence the current GRB sample. The 
well-known "Amati" correlation for long-duration GRBs be- 
tween isotropic-equivalent radiated energy Eiso, describing 
the intensity of the burst, and the photon energy at which 
the time-averaged spectrum peaks i?p,i, although proving to 
be quite robust (Amati et al. 2008; Amati 2010), has shown 
evidence of variation with redshift (Li 2007) and suscepti- 
bility to detector threshold selection effects (Butler et al. 
2007). Petrosian et al. (2009) find evidence for evolution of 
the GRB peak luminosities, but this should not affect the 
Hubble diagram, since it is the luminosity relations which 
should give the right distances for placement on the Hubble 
diagram. 

It can be argued that the kind of relativistic and geo- 
metric effects that underlie the luminosity relations should 
not be greatly affected by evolution or the metallicity of the 
progenitor (Schaefer 2007). It is conceivable that a better 
understanding of GRB physics in the future will allow them 
to be used as "standardizable" candles, and put their utility 
for cosmological parameter estimation and discrimination 
between cosmological models on a firmer basis. Ongoing ob- 
servational programmes such as Swift continue to contribute 
to this aim. We need to know more about the physics of 
the GRBs, and we need more high-quality measurements of 
GRB redshifts, light curves and spectra. 

In the meantime, however, we can obtain glimpses of 
the potential applications of standard candles whose range 
extends into the era of decelerating cosmic expansion. In 
the present study, the correlations are forced to be a good 
fit by incorporating the additional "systematic error" term, 
computed such that it makes the x^ of the best-fit correla- 
tion equal to one. This term contributes (in quadrature) to 



the uncertainty in the log of the isotropic luminosity which 
propagates through to the Hubble diagram x^ via eqs (12) 
and (13). A single GRB at a redshift of 5 or 6, with better- 
determined physical characteristics, potentially carries more 
statistical power than a single SN la at z = 1.7 because of 
the Hubble diagram "lever arm" — the Hubble diagram at 
redshifts z > 2 changes with a different cosmological model 
or cosmological parameters much more than it does at lower 
redshifts. We obtain results that are not inconsistent with 
current models, with certain acknowledged caveats. In par- 
ticular, there is much scope for progress in improving the 
GRB Hubble diagram, and much to be gained. 
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APPENDIX A: TIMESCAPE CALIBRATION 
CURVE EQUATIONS 

The five calibrations in figures 2-6 are based on the bisector 
of the two ordinary least squares fits (Isobe et al. 1990). For 
the ith luminosity indicator, the best-fit line has the form 
Yi = a + bXi, where Xi — log(indicator) ± log(l -I- z), where 
the sign of the redshift factor depends on the indicator. The 
five best-fit lines for the TS model (those for which the HD 
X^ is a minimum), their associated uncertainties, and the 
corresponding ACDM values for a and b, are given below. 



(i) Lag time vs. Luminosity: 
For the timescape model: 
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(ii) Variability vs. Lum,inosity: 
For the timescape model: 
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log L = 52.57 + 1.80 log 



\V{1 



L 0.02 



<^lo 



+ 



/0.43436o-v\2 2 

' ' +f^V,sys, 



(A3) 



(A4) 



V 






where aa ~ 0.34, at — 0.20, and ay ^y^ = 0.35 gives a 
reduced x^ of one. For the ACDM calibration, we find a — 



52.50 ± 0.34, b = 1.77 ± 0.20 and a^^^^^ 



0.35. 
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(iii) -Epeak vs. Luminosity: 
For the timescape model: 

log L^ 52.31 + 1.71 log [^.^^^ 



O"log L — 



O-a + jo-ftlog [- 



^pcak 



(1 + ^) 



300 keV 



ir 



(A5) 
(A6) 



fOAUSbaE^^2 ^ 

where aa = 0.24, ctj, = 0.10, and cr|; ^.^j^.sys = 0.34 gives 
a reduced x^ of one. For the ACDM calibration, we find 
a = 52.24 ± 0.24, b = 1.69 ± 0.10 and o-l^^^^.y, = 0.34. 



(iv) Speak VS. E-,: 

For the timescape model: 

log E^ = 50.64 + 1.63 log [ 



f^loK E., - 



'., = 0-a + jo-fclog 1^ 



-Epoak(l + Z) 

300 keV 

gpcak(l + ^)ll^ 



^]} 



+ 



0.4343&aBp„,,^^ , _2 



S, 



peak 



300 keV 

2 

+ ""£;_,, sys I 



(A7) 



(A8) 



where aa = 0.28, Ub = 0.10, and a% ^y^ = 0.17 gives a 
reduced x^ of one. For the ACDM calibration, we find a = 
50.58 ± 0.28, 6 = 1.62 ± 0.10 and cr| , ,y, = 0.15. 



(v) Rise time vs. Lum,inosity: 
For the timescape model: 



log L = 52.61 - 1.25 log riHiU 



'rRT(l + zY 
s 



""log L 



"» +h^°gi 0.1 s J I 



(A9) 
(AlO) 



0.43436o-Ep^^, 

-C/pcak 

where Ga ~ 0.11, ah = 0.11, and cr^j^^^sys = 0.48 gives a 
reduced x^ of one. For the ACDM calibration, we find a = 
52.54 ± 0.11, b = -1.23 ± 0.11 and cr^j^^-sys = 0-47- 
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